srs.ratio2=function(data,x,y,n,mux){
  data=as.name(data); x=x; y=y; n=n; mux=mux
  plot(x,y,main=data)
  r=sum(y)/sum(x); s2r=sum((y-r*x)^2)/(n-1)
  vhat.r=(1/mux^2)*(s2r/n)
  bound=2*sqrt(vhat.r)
  lower=r-bound; upper=r+bound
  cat("","\n","Results from SRS: Data =",data,"\n",'Estimation method = Ratio with unknown N','\n',"n =",n,"\n","Ratio r =",r,"\n","Vhat r =",vhat.r,"\n","Bound =",bound,"\n","Lower Bound =",lower,"Upper Bound =",upper,"\n","")
  results=list(r=r,data=data,n=n,vhat.r=vhat.r,bound=bound,lower=lower,upper=upper)
}
# to use the function with its call:
# srs.ratio2(data,x,y,n,mux)

# data: name of dataset, in quotes
# x: vector (data)
# y: vector (data)
# n: sample size
# mux: value of mean of X